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Abstract 

The hard square model in statistical mechanics has been investigated for 
the case when the activity z is —1. For cyclic boundary conditions, the 
characteristic polynomial of the transfer matrix has an intriguingly simple 
structure, all the eigenvalues x being zero, roots of unity, or solutions of 
x 3 = 4 cos 2 (irm/N) . Here we tabulate the results for lattices of up to 12 
columns with cyclic or free boundary conditions and the two obvious orien- 
tations. We remark that they are all unexpectedly simple and that for the 
rotated lattice with free or fixed boundary conditions there are obvious likely 
generalizations to any lattice size. 
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1 Introduction 

The hard square lattice model has long been the subject of study. In 1966 
Runnels and Coombs pQ performed numerical calculations on lattices of up 
to 24 columns (an amazing achievement for the computing power then avail- 
able). Gaunt and Fisher had previously obtained high and low-density se- 
ries expansions in 1965 [2J, and these were extended by Baxter et al [3J. For 
the case z = +1, Baxter also calculated the free energy to 43 digits of accu- 
racy using the "corner transfer matrix" technique. [4 J Recently Fernandes et al 
have made numerical studies of various hard core two- and three-dimensional 
systems. [6] 

All these calculations naturally focussed on cases of physical interest, when 
the activity z is real and positive, and were particularly concerned with lo- 
cating and investigating the critical point at z ~ 3.796, where the system 
undergoes a transition from a fluid to a solid state. However, as a result of 
their analysis of the series expansions, Gaunt and Fisher did predict singular- 
ities on the negative real axis, at z — —0.12 and z ~ —3.80. Between those 
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values the free energy may have branch cuts or further singularities, and may 
well be sensitive to the boundary conditions imposed. 

Fendley et al [7] studied the case when z = — 1, using the usual cyclic 
boundary conditions and orientation. For lattices of N columns, with N < 15, 
they found that the characteristic polynomial Pn( x ) = det(xl — T) of the row- 
to-row transfer matrix is remarkable simple, being a product of factors of the 
form x m db 1. This implies that all the eigenvalues of T lie on the unit circle. 
They conjectured that this behaviour held for lattices of any width. 

This conjecture was verified by Jonsson,[8] who established an equivalence 
between the model and the counting of particular rhombus tilings of the plane. 
Moreover, he obtained an algorithm for counting these tilings that grows only 
polynomially (instead of exponentially) with the size of the lattice. He was 
therefore able to calculate Pn{%) explicitly for lattices of up to 50 columns. 

More recently, Jonsson [9] has investigated the effect of first rotating 
the lattice through 45° before imposing toroidal boundary conditions, and 
Bousquet-Melou et al^0\ have looked at the problem using free boundary 
conditions. Again they found that Pn{x) is a product of simple factors. 

We have explicitly calculated Pn(%) for N < 12, and for completeness 
present the results here for the three cases mentioned, plus the remaining 
case of the 45° orientation with free (or fixed) boundary conditions. This last 
case is particularly simple, leading us to conjecture its form for general N. 

2 The 90° orientation 

Fendley et al [7J considered the square lattice C drawn in the obvious way as 
in Figure [H with M rows and N columns and R = MN sites. For the hard 
squares model, with each site i is associated an occupation number <jj, which 
takes the values or 1, corresponding to whether the site is empty or contains 
a particle. At most one particle can lie on any site. 
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Fi gure 1: The square lattice of solid lines and circles, showing two adjacent sites i , j. 
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Particles are also not allowed to occupy adjacent sites. We can express 
this by defining a weight function W(ai,aj) for each edge of C by 

W(a,b) = 1 if a,b not both 1 , (2.1) 
= if a = b = 1 . (2.2) 

Let <r= {<7i, . . . , <tr} denote the set of the occupation numbers of all the 
sites of £, and let 

n(<7)=£>; , (2.3) 
the sum being over all the R sites i. Then the partition function is 

Zm,n = 5> n(<T) UWfaaj) , (2.4) 
°" <<j) 

the product being over all edges of £ and the sum over all values of 

(71,(72, • • • )<7K- 

Define the partition function per site 

K = K (z) = (Z m ,n) 1/R ■ (2.5) 

The parameter z is known as the activity. For real positive values of z we 
expect k(z) to tend to a positive limit as M, N — > oo, to be independent of 
the order in which this occurs, and indeed of the boundary conditions imposed. 
The free energy or entropy is proportional to — \og k(z) and is the function 
investigated by Runnels and Coombs, p] Gaunt and Fisher, [2] and Baxter et 

al pg. 

If we impose toroidal boundary conditions, so that row M is followed by 
row 1, and column N by column 1, then it is straightforward to show that 

Z m ,n = trace (T N ) M , (2.6) 

where T/v is the transfer matrix, with entries 

N 

{T N )*,<r> = n^%^i+i)%^i) (2-7) 

J=l 

where now a = {a\, . . . , 0^}, a' = {a[, . . . , a' N } and crjv + i = o\. The RHS 
is the combined weight of two adjacent rows of C, with occupation numbers 
o"i, . . . , ctjv in one row, and a[, . . . , a' N in the row above, as in Figure [2j 

As defined here, Tjy is a 2^ by 2^ dimensional matrix. However, many of 
its rows vanish, corresponding to two horizontally adjacent occupation num- 
bers (Tj,<Tj+i being both unity. It is natural to exclude such rows and to 
impose the condition 

(7j , <Tj+i not both one, j = 1, . . . ,N (2-8) 

on both the sets a and a'. This does not change the RHS of (|2.6p . If S is the 
two-by-two matrix 
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Figure 2: Two rows of £, corresponding to the element T a , a i of the transfer matrix T. 




then the matrix T/v is of dimension 

Dn = trace S N . 

The characteristic polynomial of the matrix TV is 

p N (x) = det (x I- T N ) . (2.9) 

For arbitrary values of the activity z, one would not expect this to have any 
simple structure. However, the case 

z = -l (2.10) 

seems to be very special. Fendley et al [7j evaluated Pn{x) for N = 1, . . . , 15 
and found that it to be amazingly simple, being in each case a product of 
factors of the form 

(x n ±l) . 

Thus all the eigenvalues of T/v lie on the unit circle. 
We define the functions 

fn = fn(x) = X n - 1 , (2.11) 

and give Fendley et aFs results up to N = 12 in Table [TJ In each case the 
denominator is a factor of the numerator, so each expression is in fact a 
polynomial in x, of degree Djy. 
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Table 1: Pn( x ) for the 90° orientation with cyclic boundary conditions. 



Free and fixed boundary conditions 

One does not have to impose toroidal (i.e. cyclic in both the vertical and 
horizontal directions) boundary conditions on the lattice C. Other boundary 
conditions are of interest in statistical mechanics, in particular cylindrical 
boundary conditions, where £ begins on the left at column 1 and ends on the 
right at column N. If one still imposes cyclic boundary conditions from top 
to bottom (so that row M is followed by row 1), then Zm n is again given by 
(I2.6D and (|2.7j) . except that the single factor W(on, ctv+i) is omitted. 

This is equivalent to requiring that (|2.8p hold only for j = 1, . . . ,N— 1, so 
now T/v is a square matrix of dimension 

i=i j=i 

This case is discussed by Bousquet-Melou et q/|10|. who calculate Pn(x) 
for iV = 1, . . . , 10. We have extended these calculations to N = 1, . . . , 12 
and give the results in Table [2 Again the denominators are factors of the 
numerators, so each expression is a polynomial, for instance 

Pi(x) = x 2 -x + 1 , P 2 (x) = (x - l)(x 2 + 1) . 

We observe that, surprisingly, these results are similar in form to those for 
the cyclic case, in particular all the eigenvalues, i.e. the zeros of Pn(x), lie on 
the unit circle. 

We can also consider fixed boundary conditions, for instance we could fix 
o"x, o"tv, cr[, cr' N in (|2.7p to be zero. (The factor W(aiy,ai) is then unity, so 
it is irrelevant whether it is removed.) However, since then W(ai,o~2) = 
W(ctn-i,o~n) = 1 for all 02 and ctn-i, this is merely the same as imposing 
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Table 2: Pn(x) for the 90° orientation with free boundary conditions. 



free boundary conditions on columns 2 and N— 1 of £, so for the 90° orientation 

[PN(x)] SxBi = [P N - 2 (x)] bee . (2.12) 

3 The 45° orientation 

We now consider the square lattice C turned through 45°, as in Figure El 

If we sum over the occupation numbers on alternate rows (the solid circles 
in the figure), then we obtain a new square lattice C\, shown by the open 
circles and dotted lines in Figure O which has the usual 90° orientation used 
in the previous section. The four sites round a face of C\, e.g. those with 
occupation numbers a, b, c, d shown, have the weight 

1 

Wi(a,b,c,d) = ^z e TT(a,e)W(6,e)lT(c,e)W(d,e) . (3.1) 

e=0 

Thus 

Wx(a,b,c,d) = l + z ifa = b = c = d = , (3.2) 

= 1 otherwise (3-3) 

and we see that z = — 1 is special in the sense that Wi(0, 0, 0, 0) then vanishes. 
The definition (|2.4p is equivalent to 

Z m ,n = $> n(<7) I[W 1 (a i ,a j ,a k ,a l ) , (3.4) 
a 

the product being over all faces (i,j,k,l) of C\ and the outer sum over all 
values of the occupation numbers a on L\. Again n(<r) is defined by (|2.3p . 
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Figure 3: The square lattice C turned through 45°, indicated by all the circles and the solid lines, 
and the reduced lattice L\ shown by open circles and dotted lines. 

but now the sum is over all sites of L\ . We now take M, N, R to be the number 
of rows, columns and sites of C\, so again R = MN. 

The transfer matrix is no longer given by (|2.7p . but by 

N 

(TN)a,a> = Y[z a3 W 1 (a j ,a j+1 ,a' j+1 ,a' j ) (3.5) 

3=1 

and is of dimension 2^. None of the rows or columns now vanish, but many 
of the rows of Tjv are equal or opposite. This is because two rows with 
a = {a, . . . , b, 1, 0, 1, c, . . . , d} and {a, ... ,b, 1,1,1, c, ... , d} must be equal and 
opposite. (The central occupation number is different, but the face weights 
on either side are 1 for each of the two cases.) Thus many of the eigenvalues 
of Tat will be zero, and we do observe this. 
For positive integers n, define 

9n = 9n{x) = x 3 - n . (3.6) 

If we impose cyclic boundary conditions, so that column iV of C\ is followed by 
column 1, then for z = — 1 and N = 1, . . . , 12 we find that the characteristic 
polynomial P/v(x), as defined by (|2.9p . is given as in Table [3j 

If N is divisible by 3, we see that there are extra factors whose zeros do not 
in general lie on the unit circle. Nevertheless they do have a simple structure: 
in each case Pn{x) is a product (or ratio) of factors x, f n and of 

s(m,N) = x 3 - 4cos 2 (vrm/iV) . (3.7) 

To see this, note that g 4 = s(0,N), g 3 = s(N/6,N), g 2 = s(N/A,N), g x = 
s(N/3,N) , and 

x 9 - 6x 6 + 9x 3 - 1 = s(l, 9) s(2, 9) s(4, 9) , 
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Table 3: Pn( x ) for the 45° orientation with cyclic boundary conditions. 



x 6 -4x 3 + 1 = 8(1,12)8(5,12) . 

This is the case investigated by Jonsson in his second paper [9j. Using the 
equivalence to rhombus tilings of the plane, he does indeed show that such 
factors s(m, N) occur when N is divisible by three. 

Thus the eigenvalues of T/y lie either at the origin, on the unit circle, or 
their cubes lie on the interval (0,4] of the positive real axis. 



Free and fixed boundary conditions 

As with the 90° case, we can instead impose free boundary conditions, so the 
factor Wi((tn, 01, a'±, a' N ) (occurring when j = N) in the product in (|3.5p is 
omitted. Now we find, for N = 1, . . . , 12 that Pn(x) has the very simple forms 
given in Tabled! 

For a given N, define m by 

m = (4 fc -l)/3 ifiV = 3A; (3.8) 
if N = 3k + 1 (3.9) 

= (l + 2 2fc+1 )/3 ifiV = 3A; + 2 , (3.10) 

where k is an integer. Also set 

r = 1 — mod (N, 3) , n = 2 N — 3m — r 

(so r = 1,0 or —1). Then the results of Table H] are fitted by 

p N ( X ) = ip N ( X )] hee = X "f™r 1 . (3.11) 

As with the 90° orientation, we can also consider fixed boundary conditions, 
in particular by fixing the occupation numbers in the first and last columns 
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Table 4: Pn(x) for the 45° orientation with free boundary conditions. 



of C\ to be zero. However, if we go back to the original lattice C (the solid 
lines in Figure [3]) which has 2N — 1 columns, this is equivalent to imposing 
free boundary conditions on columns 2 and 2N — 2. This means that 

[ZmMbxbA = l Z M,N-i] {ree ■ (3.12) 

The same relation does not quite hold for the transfer matrices, because XJv for 
fixed boundary conditions is of dimensions d, while T/v_i for free boundaries 
is of dimensions 2d, where 

d = 2 N ~ 2 . (3.13) 

For fixed or free boundary conditions we see from Figure [3] that there are 
two sorts of rows of the original lattice C: the ones shown as open circles 
containing N sites per row, and the ones shown as solid circles containing 
only N — l sites. This ensures that the transfer matrices factor: 

[TjV-l]f r ee = $nRn , [T/v] fixed = RnSn , (3-14) 

where R n is a d by 2d matrix and Sn is 2d by d. This means that d of the 
eigenvalues of T/v_i(free) vanish trivially, while the other d are the same as 
those of T/v (fixed). Hence 

[PN(x)] &xcd = [PiY-i(x)] free /x d . (3.15) 

Thus for both orientations the characteristic polynomial for fixed boundary 
conditions is obtainable from that for free boundary conditions. 

4 Conclusions 

For cyclic boundary conditions, with either the 45° or 90° orientations, Jons- 
son has shown that Pn{x) is product or ratio of factors of the form x, x m — 1, 
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and (for the 45° case with iV divisible by three) x 3 — 4cos 2 (7rm/A r ). Some- 
what remarkably, these properties appear to remain true when we impose 
free or fixed boundary conditions. In particular, we conjecture that for the 
45° free or fixed cases, Pn(x) is a product or ratio only of the simple factors 
x, fi(x), fs{x), being given by (|3.1ip and (|3.15H . 

For the 90° cyclic case, if we take the limit when M, N become large while 
remaining co-prime, then Zm,n = 1- On the other hand, from our conjecture 
(3), for the 45° free case, if N = 3k + 1, then r,m in (|3.1ip are both zero, 
giving 

Z m ,n = ■ (4.1) 

From (12, 5j) k for these two cases is 1 and 0, so its limit is certainly 
dependent on the boundary conditions. This can happen, because when z = 
— 1 the Boltzmann weights are not all positive. 

So z = — 1 is a very special non-physical case of the hard squares model. 
Even so, the transfer matrix does then appear to have very intriguing prop- 
erties. 
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